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ABSTRACT 

Context. V838 Mon erupted at the beginning of 2002. In the course of the outburst the object evolved to low effective temperatures 
and declined as a very late M-type supergiant. Among various scenarios proposed to explain the nature of the outburst, the most 
promising is a stellar merger event. 

Aims. We aim at studying the structure and evolution of the object in the decline from the 2002 eruption. 

Methods. We obtained spectroscopic observations of V838 Mon in January-March 2009 with UVES/VLT. The results are analysed 
and compared with similar observations obtained in October 2005 with HIRES/Keck. 

Results. The most striking difference between 2009 and 2005 is a complete absence of the B3 V component and of the [Fell] 
emission lines in 2009. The present spectrum displays only the spectrum of the 2002 eruption remnant. It resembles that of an ~M6 
giant, although the molecular bands in V838 Mon are deeper than those in standard stellar spectra of a similar spectral class. Several 
atomic lines, which displayed P-Cyg profiles in 2005, are now dominated by pure absorptions. Some of these lines, however, show 
a narrow emission component, which, as we argue, measures the radial velocity of V838 Mon. The resulting heliocentric velocity is 
71 km s"' , which very well agrees with the velocity obtained from SiO maser observations. The atomic lines and the molecular bands 
show very complex kinematics. In some atomic lines and high-excitation molecular bands we observe matter infalling in the object 
atmosphere. The infall components were already observed in 2005, but were less pronounced and present in fewer lines than in 2009. 
We argue that the most negative radial velocity components seen in the resonance atomic lines and in the low-excitation molecular 
bands were formed in the ejecta of the 2002 eruption. The B3 V companion most probably became engulfed in an opaque dusty matter 
of the 2002 V838 Mon ejecta. 

Key words, line: identification - line: profiles - stars: individual: V838 Mon - stars: peculiar - stars: late-type 



1. Introduction (lKulkamiet al.ll2007h and NGC300OT2008 (iBond et al.ll2009l: 

fT^ 'Ber ger et al.l 20091) . probably also belong to this class of stellar 

Q The erupt ion of V838 Mon was observed in January-March eruptions 
^-H 2002 (e.g. Muna ri et al.ll2002al) . The object attracted significant 

^ ' public attention, m ainly because of its spectacular light echo The nature of the V838 Mon- type eruptions has not yet been 

L; ! (IBond etalJl2003l) . The eruption itself was also very interest- Mly explained. As discussed in [Tylenda & Soker (2006), ther- 

. ^ ' ing and intriguing. The light curve could have been taken for monuclear mechanisms such as classical nova (Iben & Tutukov| 

V that of a slow nova, but the colour and spectroscopic evolution ^1223 or a late He-shell flash (Lawlor 2005 ) should be excluded. 

Cl of the object suggested a different nature. In the course of the Instead, Tylenda & Sokeij (|2006|) showed that all the main obser- 

P3 eruption the spectrum evolved to late spectral types and the ob- vational characteristics of the V838 Mon-type eruptions can be 

" ' ' ject declined as a very late M-type supergiant. At present, soon consistently understood as resulting from stellar collisions and 

a decade after the eruption, the remnant still shows a spectrum, mergers. This idea obtained a strong support from archive photo- 

which can be classified as an ~M6 giant or supergiant. metric observ ations of the progenitor of VI 309 Sco analysed in 

A few other stellar eruptions show characteristics very sim- iTylenda et al. | (|20Tl). The data evidently show that the progen- 

ilar to those of V838 Mon. They are named "V838 Mon-type '^r was a contact binary quickly evolving towards it^erger. 

objects" or "red n ovae". Apar t from V838 Mon, t his class in- According to| Kashi et al.| (l2010|) and| Kashi & Sokei| ^1^, cer- 

cludes M31 RV (Mould et al. 1990), V4332 Sgr ( Martini et al J ^ain optical transients may result from mass transfer events m 

[1999) . and VI 309 Sco (Mason et al..,2010.) . Some extragalac- ^^"^y systems. 

tic intermediate-lumnosity optical transients, hke M85 OT2006 A B3 V companion of V838 Mon was disc overed when 

the main object declined from the 2002 eruption ( Munari et al.l 



Send ojfprint requests to: R. Tylenda |2002b). The companion remained constant in brigh tness until an 

* Based on observations made with ESO Telescopes at the Paranal eclips e-hke event in November-December 2006 (iMunari et al.1 
Observatory under programme ID 382.D-0152(C). l2007h . In 2004, narrow emission lines, predominantly from 



2 



Tylenda et al.: V838 Mon spectram in 2009 



[Fell], appeared in the spectrum of V838 Mon (iBarsukova et alJ 
12006 ^. The lines strengthened with time and reached a maximum 
near the 2006 eclipse-Hke event (Munari et al. 2007). The emis- 
sion lines were most likely excited by the radiation of the B3 V 
companion in the matter ejected by V838 Mo n in 2002 that ap- 
proached the companion (iTvlenda et al.ll200"9l) . The eclipse-like 
fading of the companion in 2006 was accordingly caused by a 
cloud of dusty matte r crossing the line of sight of the companion 
(ITvlenda et al.ll2009l) . 

A few months after the eclipse-like event, the B3 V com- 
panion started a new decline, this time slower but significantly 
deeper and more complexQln contrast to the 2006 event, the sec- 
ond decline of the companion was also accompanied by a decline 
of the [Fell] emission lines. This indicated that the radiation of 
the B3 V companion faded not only for the observer but also for 
the matter flowing by the companion. Thus the event cannot be 
explained by another dusty cloud crossing the line of sight. The 
companion was instead engulfed in dense and opaque matter 
Two years after the beginning of the decline, last traces of the 
B3 V companion disappeared from the spectrum of V838 Mon. 

In March 2008, a field of V838 Mon was observed by 
XMM/Newton and a varia ble X-ray source w as detected near 
the position of the object (Anton ini et al.l20Toh . The X-rays may 
have resulted either from magnetic activity of the merger rem- 
nant or from interactions of the 2002 ejecta with the B3 V com- 
panion. No such source was detected in Chandra observations 
obtained a year after the outburst, nor in January 2010. 

In October 2005, we obtained a high-resolution spectrum 
of V8 38 Mon with the HIRES/Keck instrument ( Kamihski et ak 
|2009(). A detailed analysis of that data was made in Tvlen da et aP 
(|2009|) . In the present paper, we report on similar observations 
fliat were obtained in January-March 2009 with the UVESA'LT 
instrument. 



2. Observations 

2.1. Spectroscopy 

V838 Mon was observed with the UVES spectrograph 
jDekker et al.l2000l) on the Very Large Telescope (VLT) in 2009 
in three observing runs: on 2 January, 13-14 February, and 22- 
23 March. For each run two exposures were obtained with a to- 
tal exposure time of 85 min. Technical details of the individual 
exposures, i.e. observation date, UT time and Julian date at the 
exposure start, integration time, averaged airmass, and average 
seeing conditions calculated over the exposure time, are given in 
Table[T] The observations were performed in the UVES standard 
red arm 580 setting, which gives a spectral coverage of 4770- 
6825 A with an inter-chip gap between 5785 A and 5820 A. 
We used a slit with a width of 0.7 arcsec and the UVES im- 
age slicer no. 1. These settings provided a spectral resolution of 
about A/AA^52 000. 

All spectra were reduced and cal ibrated in IRAPEI using stan - 
dard procedures for Echelle spectra (' Willmarth & Barneslll994l) . 
The reduction procedures included bias- and flat-field correc- 
tions, removal of scattered light, and spectra extraction with 
background subtraction. Cosmic ray events were rejected from 
the CCD frames with the Laplacian cosmic ray identification 

' see http://jet.sao.ru/jet/~goray7l 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



algorithm (Ivan Dokkumll200TI) . Observations were wavelength- 
calibrated using spectra of Th-Ar lamps, and dispersion solu- 
tions were fitted to identified lines with a typical accuracy of 
rms=3 mA. Flux calibration was based on observations of the 
spectrophotometric standard star Hiltner 600 and its spectropho- 
tometric fluxes tabulated every 50 A in Ha muv et al.l (1 19941 
il992i) . We asses the overall accuracy of the relative flux cali- 
bration to be about 5%, but it is slightly poorer within strong 
telluric absorption bands and in spectral regions correspond- 
ing to the broad B aimer lines in the spectrum of the standard 
star. The flux calibration of the spectra from 24-25 March was 
based on standard-star observations on 26 March and therefore 
the absolute flux calibration of the March spectra is most uncer- 
tain. All spectra were shifted to the heliocentric rest frame, and 
all Echelle orders from two expositions obtained in the same 
run were merged into one spectrum. The signal-to-noise ratio 
reaches a value of 120 in the spectral region with the highest 
flux level. 

In order to check the accuracy of the absolute flux calibration 
of our spectra, we estimated their corresponding V magnitudes 
by integrating the spectra convolved with the response curve of 
the V filter The inter-chip gap was interpolated w ith the shape of 
the spectrum reported in iKaminski et al.l (|2009) as a reference. 
The magnitudes derived from the spectra obtained in January, 
February, and March are 15.4, 15.5, and 15.3, respectively. They 
can be compared to the V values of 15.36 and 15.34 derived 
from our photometric measurements made in December 2008 
and March 2009, which we described in Sect. l2.2l (see Table|2). 

The three spectra from January, February, and March have 
essentially the same shape. However, they are systematically 
shifted towards each other on the flux scale. The flux level in 
the spectrum from February is on average 7+2% lower than that 
of the spectrum from January, and the last spectrum from March 
is 14+2% brighter than the spectrum from January. These flux 
variations are within the accuracy of the absolute flux calibration 
of our spectra. We can conclude that the object had essentially 
the same flux level in the three observations. Minor differences 
between the three spectra related to individual features are dis- 
cussed in Sect. 13.4! We have summed up the three spectra and 
presente and analyse the resulting spectrum. 

2.2. Photometry 

V838 Mon was measured in the optical photometric bands sev- 
eral times in 2008-2009. On 18 December 2008 the object was 
observed in the V and Ic filters using a CCD camera at the 
1 m telescope of the South African Astronomical Observatory 
(SAAO). Photometric observations were also carried out on 
16 March 2009 in the BVRcIc filters with the 1 m SMARTS tele- 
scope of the Cerro Tololo Inter-American Observatory (CTIO). 
Finally, a BVIc photometry was performed on 5-6 October 2009 
with the 1 m SAAO telescope. 

All the above observations were reduced in IRAF using stan- 
dard procedures. The magnitudes of V838 Mon w e re deri ved 
using photometry of field stars from iMunari et al ] ('2002a'). It 
should be noted, however, that t he quality of the ph otometry 
in Munari et al. was questioned in lAfsar & Bondl (l2007l) . which 
could introduce some extra systematic errors to our photometry. 
Magnitudes corresponding to different exposures from the same 
observing night were averaged for each filter The results are 
given in Table |2] where errors are also provided; they include 
the 3xrms (root mean square) scatter in magnitudes of differ- 
ent exposures during the night and the 3xrms uncertainty in the 
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date 


UT 


JD-24e5 


object 


exp. time (s) 


airmass 


seeing (") 


2009-01-02 


05:04:29 


54833.211 


V838 Mon 


2105 


1.07 


1.53 


2009-01-02 


05:56:33 


54833.248 


V838 Mon 


3005 


1.10 


1.57 


2009-01-02 


05:45:50 


54833.240 


Hilt 600 


305 


1.17 


1.52 


2009-02-13 


01:29:13 


54875.062 


V838 Mon 


3005 


1.09 


1.57 


2009-02-14 


03:54:36 


54876.163 


V838 Mon 


2105 


1.19 


1.50 


2009-02-16 


02:58:02 


54878.124 


Hilt 600 


305 


1.18 


0.97 


2009-03-22 


01:32:37 


54912.064 


V838 Mon 


3005 


1.19 


0.86 


2009-03-23 


01:25:30 


54913.059 


V838 Mon 


2105 


1.18 


0.63 


2009-03-23 


02:06:34 


54913.088 


Hilt 600 


305 


1.48 


0.74 



Table 2. Results of the photometric measurements of V838 Mon 



in 2008-2009. 


Date 


JD 


Band 


Mag 


Error 


Observat. 


18 Dec. 2008 


2454819.45 


V 


15.36 


0.02 


SAAO 


18 Dec. 2008 


2454819.45 


k 


10.08 


0.10 


SAAO 


13 Feb. 2009 


2454876.40 


J 


7.43 


0.13 


TCS 


13 Feb. 2009 


2454876.40 


H 


6.25 


0.24 


TCS 


13 Feb. 2009 


2454876.40 


Ks 


5.50 


0.07 


TCS 


16 Mar. 2009 


2454906.52 


B 


19.11 


0.13 


CTIO 


16 Mar. 2009 


2454906.53 


V 


15.34 


0.05 


CTIO 


16 Mar 2009 


2454906.54 


Rc 


12.63 


0.06 


CTIO 


16 Mar 2009 


2454906.55 


Ic 


10.28 


0.12 


CTIO 


05 Oct. 2009 


2455109.64 


B 


18.58 


0.66 


SAAO 


05 Oct. 2009 


2455109.63 


V 


15.74 


0.28 


SAAO 


05 Oct. 2009 


2455109.63 


Ic 


10.02 


0.15 


SAAO 


06 Oct. 2009 


2455110.63 


B 


18.69 


0.43 


SAAO 


06 Oct. 2009 


2455110.62 


V 


15.65 


0.39 


SAAO 



transformation factors. In 2008 V838 Mon became the brightest 
Ic source in the field of Munari et al., making its calibrations 
uncertain in this band. 

Near-infrared (NIR) photometric observations of V838 Mon 
in the JHKs bands were caiiied out on 13 February 2009 with 
the CAIN-3 camera installed on the 1.52 m Carlos Sanchez 
Telescope (TCSfl The date of the NIR observations is the same 
as that of the second spectroscopic observing run at UVES/VLT. 
Data were reduced using CAIN data reduction script^ devel- 
oped in IRAF by the Instituto de Astrofisica de Canaiias. The 
instrumental magnitudes were converted to the standard system 
using observations of two standard stars FS 11 and FS 18. Their 
NIR magni tudes were taken from a hs t of the UKIRT faint stan- 
dard stars (ICasah & Hawardenlll992l) . The resulting averaged 
NIR magnitudes for V838 Mon are presented in Table |2] 

3. Results, analysis, and comparison to the 
spectrum obtained in October 2005 

3.1. Spectral features 

A high-resolution spec trum of V838 Mon in its post-outburst 
phase was presented in iKamiiiski et a l.' ('2009') and analysed in 
[Tylenda et al. (200§. This was a HIRES/Keck spectrum ob- 
tained in October 2005, i.e. a year before the eclipse-li ke event 
observed in November/December 2006 (see, e.g. Munari et al.l 
l2007h . A de tailed atlas of spectral features in that spectrum can 
be found in iKamiiiski et alj ( l2009h . Most of the identified fea- 

^ The 1.52m Carlos Sanchez Telescope is operated on the island 
of Tenerife by the Instituto de Astrofisica de Canarias in the Spanish 
Observatorio del Teide. 

http://www.iac.es/telescopes/cain/reduc/caindrhtml 



tures there are also present in the 2009 UVES/VLT spectrum 
analysed in this paper The 2009 spectrum is compared to the 
2005 one in Fig\t\ 

The complete absence of the spectrum of the B3 V compan- 
ion and of the [Fe II] emission lines in the 2009 spectrum is the 
most striking difference between the two spectra. These spec- 
tral components were best observed in the blue part of the 2005 
spectrum, which, i.e. A < 4750 A, was not observed in 2009. 
Nevertheless, from the available data we can quite firmly con- 
clude that the B3 V component disappeared. Firstly, in 2009, 
we do not see any trace of the broad H/3 absorption line, which 
was prominent in 2005 (see Fig[T]i. Secondly, in the deep, heav- 
ily saturated TiO bands, i.e. at -4950 A, 5850-5910 A or 6160- 
6230 A, we saw a residual flux in 2005, which was nicely in- 
tepreted as the continuum from the B3 V companion. In the 2009 
spectrum the flux in these regions goes practically to a zero level, 
i.e. significantly below the 2005 level (see Fig. [TJ. Comparing 
the observed fluxes in these regions, we can conclude that in 
2009 the B3 V spectral component was at least 30 times fainter 
than in 2005. 

As mentioned above, the emission features, which were 
present in 2005, are not seen in the 2009 spectrum. We have 
not found any remnant emissions at the positions of the [Fe II] 
lines, which were strong in 2005. Comparing the observed fluxes 
in the regions of the strongest [Fell] emission lines in 2005, we 
can conclude that these lines in 2009 are at least 20 times fainter 
than in 2005. 

The emission components of the resonance lines showing 
P-Cyg profiles in 2005 also became very faint or disappeared 
completely in 2009. Similarly, the emission band of ScO near 
6460 A is not seen now. The only emission feature, which is 
clearly present in the 2009 spectrum, is that of the Ha line (see 
Figs.|I]and[3}J). 

Almost all the atomic absorptions found in 2005 are present 
in the 2009 spectrum. Their profiles, however, changed consid- 
erably (see Sect. 13.3b . In addition, the present spectrum displays 
several absorption features of Fe I, Ti I, V I, and CrI, which were 
very weak or absent in 2005. These lines belong to multiplets 
with lower energy levels of ^1.1 eV. 

In 2009, V838 Mon shows a very rich molecular spectrum 
in ab sorption. Practically, a ll the features identified in the at- 
las of lKamiiiski et alj j2009l) are present in the UVES spectrum. 
Some of the bands are more clearly seen in the 2009 spectrum, 
which allowed us to improve their identification. The shapes and 
strengths of some of the bands have changed and this is de- 
scribed in detail in Sect. 13.61 In general, many of the strongly 
saturated bands of TiO are now seen to be less saturated. 

Certain spectral regions of the UVES spectrum were not cov- 
ered in the 2005 spectrum. Identifications of spectral features in 
these regions can be found in Fig. [1] These ranges give access 



Tylenda et al.: V838 Mon spectram in 2009 



0.8 
0.6 
0.4 
0.2 


0.8 
0.6 
0.4 
0.2 




.T 


iO a(2,0) 
[Fell] 


1 1 1 1 

TiOa(3,l) 


1 1 1 1 

AlO B-X(0,0) 


1 1 1 1 1 1 

A10B-X(2,2)aiob-X(3,3) 


Bal 


' 1 

CrI 


,TiOa(l,0) 
1 [Fell] 


' 1 ' 

TiOa(2,l) Ft 
[Fell] 
Mull "^'^ 


jj TiOa(3,!2> 
[Fell] TilJ^ 






^ [Fell] 


TiO a(4°2?J'~"''**-' 




TiO a(5,3) 




Fel 


Til 


Til Tif 


Sin+[FeII] 1 ■ 






4800 


4850 




4900 




4950 


5000 


5050 



■ Til CrI 



A10B-XAv=-l Til 



Sill . Fel . [Fell] 




VO C-X(2,0) P?' [Fell] VO C-X Cri 

Crl+Fel JlrtlCrll '■■ 

TiOf-a 




6600 



6650 



6700 
Wavelength (A) 



6750 



6800 



6850 



Fig. 1. Average spectrum obtained with UVES/VLT in January-March 2009 (black Une) is compared to the spectrum obtained with 
HIRES/Keck (grey line) in October 2005. The spectra are deredenned with Eb-v - 0.9 and Ry - 3.1, and smoothed. A synthetic 
spectrum of a B3 V star, which contributes to the 2005 fluxes, is shown for comparison (blue line). The UVES spectrum was shifted 
in fluxes by the values corresponding to the average contribution of the B3 V companion in the 2005 spectrum; the value of the shift 
is 0.21 for the two top panels, and 0.10 for the other panels. Identified atomic spectral features are indicated by red markers, while 
molecular band-heads are assigned green markers. The ordinate units are lO^'^* erg s"' cm^^A"'. 



to two strong lines of BaII/1/14934 and 6496, and absorption 
bands of AlO (the B-X system) and TiO (the y' system). 

3.2. Spectral energy distribution 

We compared the 2009 spectr um with UVES spect ra of red gi- 
ants from the Paranal Atlas (iBagnulo et al. I l200l . When the 
whole spectrum is taken into account, the best fitting spectra 
have types of M6+2. Although the reference spectra reproduce 
the spectrum of V838 Mon quite well in the regions with a 
strong pseudo-continuum, numerous absorption bands are sig- 



nificantly deeper in V 838 Mon. The same was observed in 2005 
dTvlenda et alj|2009h . 

As can be seen from Fig. \T\ V838 Mon in 2009 was sig- 
nificantly brighter than in 2005. After removing the flux of the 
B3 V star from the 2005 spectrum, we found that the flux of 
V838 Mon itself in the observed spectral region increased by a 
factor of ~2 between 2005 and 2009. This growth in the optical 
flux was caused by slight increase of V838 Mon in the effec- 
tive temperature and luminosity (see below), as well as a lower 
blanketing in the molecular bands in 2009 (see Sect. 13.6b . 
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Fig. 2. Standard supergiant photometric spectrum (M6.3 red- 
dened with Eb-v - 0.9: full line) fitted to the photometric re- 
sults obtained in December 2008 - March 2009. Full points: our 
data presented in Table |2] asterisks: mean values from the data 
compiled by V. Goranskij. See text for more details. 



We fitted standard broad-band photometric spectra to our 
photometric measurements obtained in December 2008 - 
March 2009 supplemented with mean values taken from the 
measurements obtained in the same time period that were com- 
piled by V. Go ranskij' (for details of the fitting procedure see 
lTvlendall2005h . Figure |2] shows the best fit of supergiant (lu- 
minosity class I) standards reddened with Eb-v - 0.9 to the 
photometric data. The spectral type of the fitted supergiant is 
M6.3, which agrees well with the spectral type from comparing 
the Paranal spectra to that of V838 Mon, derived above. This 
spectral type (M6.3I) corresponds to an eff'ective temperature 
of ~3270 K (Levesque et al. 2005). For a distance of 6.5 kpc 
dAfsar & Bond 2007; Kamihski et al. 2007; Sparks et al. 2008) 
the object has an eff'ective radius of ~380 Rq and a luminosity 
of ~1.5 10* Lq. Practically the same parameters are obtained 
when fitting standard giant (luminosity class III) photometric 
spectra. When com pared to the results of a similar analysis made 
in Kaminski et alJ L2009E in 2009 V838 Mon was sHghtly hot- 
ter and ~20% more luminous than in 2005. 

3.3. Profiles of atomic lines 

The atomic lines observed in our spectrum show complex struc- 
tures in their profiles. Examples are shown in Fig. |3] In the 
following, we assume that the stellar ra dial velocity is equal 
to the value derived form SiO masers (iDeguchi etal.1 l2005t 
IClaussen et al.ll2007h . i.e. Vh - 71 km s ' shown with a dashed 
line in Fig. [3] This point is further discussed in Sect. 13.51 All 



^ Note that in Kamin ski et all i200St) the effective temperature was 
de rived following the sp ectral type vs. effective temperat ure calibration 
of lSchmidt-Kaie3 ( 119821) . With the calibration of Levesq ue et al.l ( 12005!) 
the effective temprature in 2005 would have been ~3200 K. 



radial velocities given in this paper are in the heliocentric rest 
frame. 

As discussed in iKaminski et aP ( l2009l) . many resonance 
atomic lines in the spectrum of V838 Mon taken in October 2005 
showed P-Cyg profiles. All these lines (except K I /1 7698, 
Rb I /1 7800, and Mg I /1 4571, which are outside the wavelength 
range of the present spectrum) are observed in the 2009 spec- 
trum, but no more as P-Cyg ones; their emission components 
practically disappeared. The absorption component in many 
lines is now wider. It is usually more extended towards longer 
wavelengths, which is at least partly due to the disappearance of 
the emission component. Some absorptions are shifted to shorter 
wavelengths, but never beyond the maximum blueshifted ab- 
sorption observed in 2005. The maximum outflow velocity of 
^ 21 5 km s derived from the 2005 spectrum (ITvlenda et al.l 
l2009l) . remains practically the same in 2009. 

The absorption Unes show complex profiles in most cases. 
A relatively narrow absorption centred at Vh - 87 km s"' (see 
Fig. O is one of the most characteristic features seen in the 
majority of the lines. In the 2005 spect rum, this component 
was seen in the Ti I lines (IKaminski et alJl2009 ) and was inter- 
preted as a signature of an infall in the V838 Mon atmosphere 
( iTvlenda et alj|20()9l) . The wings of this feature in the Til lines 
are now observed at practically the same velocities as in 2005, 
but the feature is now deeper and the maximum absorption is at 
a slightly higher velocity. 

The infall component is the only feature seen in the VI 
lines (see panel h in Fig. [3). These are lines arising from lev- 
els with excitation energies as high as ~1.05eV. When fitted 
with a Gaussian, the component has a full-width-at-half depth 
(FWHD) of ~ 15 km s ' . We note that the V I lines were also seen 
in the 2005 spectrum, but they were then significantly weaker 
and their profil es were more noisy. T here fore, these lines wer e 
not analysed in lKamiriski et alJ (l2009h and lTvlenda et alJ (l2009l) . 
Comparing the 2005 spectrum with the present one we can now 
state that in 2005 the VI lines were of similar width but at 
slightly lower radial velocities than in 2009. This is consistent 
with the analysis of the infall component in the Ti I lines, pre- 
sented above. 

The profiles of non-resonance Fe I lines (panel e in Fig. |3]l 
show a maximum absorption at velocities consistent with the in- 
fall. However, they are wider than the profiles of the V I lines and 
considerably asymmetric. The blue wing extends well below the 
stellar velocity, indicating an outflow. Contrary to resonance Fe I 
lines (displayed in panel a in Fig. [3]) the non-resonance lines do 
not show any significant feature at Vh S 20 km s"' (structures 
seen in panel e of Fig.[3]at these velocities are caused by molec- 
ular bands). 

The infall component in the lines of Mn I and Cr I (panels 
f and g in Fig. O is separated from the blueshifted (indicating 
outflow) part of the absorption profile by a narrow emission-like 
feature peaking at the stellar velocity. In Sect. I3.5l we show that 
this is most likely a real emission feature. 

Clear outflow components, i.e. absorptions at Vh ^ 50 km s"' 
are primarily seen in lines arising from levels of low excitation 
energy, i.e. ground states (Cal, Li I, Fel: panel a in Fig.|3j Bal, 
Ball: panel b; MnLpanel f) and levels with energies $0.05 eV 
(Ti I lines in panel cjj. They are also seen in the Cr I lines (panel 
f and g), however, which arise from levels of energies as high as 



^ Th ere is a mistake in IKaminski et alj i2009h and ITvlenda et al.l 
j2009l) . where the upper energy level of the Ti I lines was given as the 
lower one. 
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Fig. 3. Sample of profiles of atomic lines in the average UVES spectrum from 2009. The vertical line marks the velocity of the 
maser of V838 Mon. Some profiles were shifted in ordinate for more clarity. Many structures, particularly seen at negative radial 
velocities, are caused by blending molecular bands. 



~1.0eV, i.e. excitation energies similar to those of the VI lines, 
which do not show any outflow (panel h). 

In the 2005 spectrum, a relatively narrow absorption com- 
po nent (NAC ) was obser ved at Vh - -82 km s"' (see Fig. 6 
in iKamiriski et alJ l2009h . It was particularly pronounced in 
Rb U 7800 and Mg U 457 1 , but also present in Mn U 5394 and 
5432. The two former lines are, unfortunately, outside the spec- 
tral range of the present spectrum, but the Mn I lines do not show 
any clear feature at the above velocity in 2009. However, an ab- 
sorption component at practically the same velocity is now ob- 
served in Ca I /1 6572 and Li I /1 6708 (see panel a in Fig. |3] and 



Fig. nil. In 2005, these two lines presented rather smooth absorp- 
tion profiles extending down to Vh - -120 km s"'. The 2009 
component is wider, however, than the similar feature in 2005, 
particularly in Li I /i 6708. 

In the 2009 spectrum, we can easily identify a very narrow 
absorption component at Vh - -92 km s"'. It is well seen in 
the Bal, Ball, and Til lines (panels b and c in Fig. [J]), but a 
trace of it seems to be also present in the blue wing of the wider 
feature in the Cal and Li I lines (panel a in Fig. |3]l. This very 
narrow feature has a Gaussian profile with FWHD - 5.6 km s 
which is the same as the FWHM of the instrumental profile. A 
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similar feature was observed in 2005 in the Ti I lines, but at Vh - 
-90 km s"'. It was also wider than now, but mostly because of 
the lower resolution of the 2005 spectrum. 

The Ca I /1 6 1 02 and 6 1 22 lines arise from the highest energy 
levels (-1.9 eVf| among the lines identified in the 2009 spec- 
trum. Consequently, they are likely to be formed deeply in the 
photosphere. This is supported by their central positions being 
close to the stellar velocity (see panel d in Fig.O. Their wings, 
extending from Vh - 10 km s"' to Vh - 100 km s"' , show an 
outflow and infall in the atmosphere of V838 Mon. 

The only clear emission line seen in the 2009 spectrum is 
that of Ha (see panel d in Fig. |3]l. The observed position of the 
line gives Vh - 29 km s The line is narrow with FWHM ^ 
25 km s"'. The approximate line equivalent width (EW) and 
dereddened flux are -0.24 A, and 3.3 x 10"'"^ erg cm^^s"', re- 
spectively. This feature was absent from the 2005 spectrum. 

3.4. Variability on a time-scale of months 

A characteristic flow time-scale of the outflow in V838 Mon 
is T ^ /?*/Vexp - 20 days (for ^ 400 R© and Vexp - 
150 km s '), which is approximately the time span between our 
three observing runs. Certain differences in individual line pro- 
files can be noticed between these three spectra. The changes are 
most pronounced in the Lil A 6707 line, which is expected to 
have the lowest optical depth among the strong absorption lines. 
The variability in Lil line is shown in the top panel of Fig. |4] 
The feature between the deep absorptions appears to slightly 
decrease with time. Interestingly, the profile of the intercom- 
bination line of Ca I /I 6572, which exhibits a similar shape as 
Lil /l 6707, has not undergone similar changes; it shows instead 
only a minor change in the March spectrum when the absorption 
profile at around the maser velocity became slightly deeper than 
in the two earlier epochs (see bottom panel in Fig. |4). Other lines 
do not show any clear signs of variability above their local noise 
levels. 

3.5. Emission in atomic lines at the stellar velocity 

The profiles of many lines displayed in Fig. |3] show an 
emission-like component at the radial velocity of the SiO maser 
(71 km s"'). This is particularly the case for the Mnl and CrI 
lines. This can be a real emission, in some sense a remnant of the 
P-Cyg emission component observed in 2005. Another possibil- 
ity is that the whole profile is purely absorptive, but the higher 
flux near the maser velocity simply means that there is little ab- 
sorbing matter at these velocities. 

The case of the Cr 1/15409 (see panel f in Fig. [3]), where 
the emission-like component peaks slightly above the continuum 
level, argues in favour of the former possibility. Additional argu- 
ments arise from comparing profiles of pairs of lines with the 
same lower energy level, but different oscillator strengths (tran- 
sition probabilities). This is the case, for instance, for Cr I /1 5346 
and /1 5349. Both lines have the same lower level, but the former 
one has a gf -value twice higher than the latter In the case of 
pure absorption, the A 5346 line is expected to be deeper than the 
/1 5349 one. Panel g in Fig. [3] shows that this is indeed the case 
except for radial velocities 550 km s '. For higher velocities, in 
particular in the vicinity of the maser velocity, the situation is 
reversed, i.e. the flux in the A 5346 line is higher than in the the 
/1 5349 one. The emission-like component is clearly stronger in 

' There is a mistake in iKaminski et alj ( I2009I) . where the upper en- 
ergy level of these Ca I lines was given as the lower one. 





0.8 






1 1 






0.7 








' yO!A ■ ■ ' ■ 


-'< 






Lil 




If \ 




0.6 












0.5 
























0.4 
























0.3 






















'o 


0.2 




\ / \. 




02Jan./1.03 ^l" 


<^ 


0.1 






13Febr. \ 








22 March/1.23 ^ 











1 
















6702 


6704 6706 


6708 


6710 6712 



Wavelength (A) 



2 0.4 



Cal 


'A A 




02Jan./1.03 




13 Febr. 




22 March/1.22 



6568 6570 6572 6574 6576 6578 

Wavelength (A) 

Fig. 4. Spectral variability in the profiles of the lAl A 6101 and 
Cal /I 6572 lines as seen in our three UVES spectra from 2009. 
The January and March spectrum were rescaled to fit the local 
pseudo-continuum in February. The feature at around 6575 A in 
the February spectrum caused by a cosmic -ray hit. The vertical 
dashed line marks the object heliocentric velocity of 7 1 km s ' . 



the former line. This shows that the emission-like component at 
Vh = 71 km s"' is indeed of emission origin. The same conclu- 
sion can be drawn when comparing the Mnl /1 5394 and /1 5432 
lines (panel f in Fig. [3]). Both lines have the ground level as their 
lower level, but the former one has a value almost three times 
higher than the latter Even the profiles of Fe I /1 5 1 10 and A 5060 
(see panel a in Fig. [3]), although they do not display any clear 
emission component, show that the high-velocity regions are af- 
fected by emission. These two lines arise from the ground state, 
but the A 51 10 has a transition probability almost 40 times higher 
than the /1 5060. For Vh % 50 km s ' the former line is consid- 
erably deeper than the latter one, but for Vh ^ 50 km s ' the 
opposite is observed. 

We can state, therefore, that an emission component is ob- 
served in some atomic lines. This is probably a residual emis- 
sion component of the P-Cyg profiles observed in 2005 and is 
most likely produced in the moving (outflowing and/or infalling) 
matter above the photosphere that is observed above the photo- 
spheric limb. The extent of the emitting region is probably small 
compared to the photospheric radius and therefore, we mainly 
observe matter moving perpendicularly to the line of sight in the 
emission component. Hence, the emission is expected to be ob- 
served at radial velocities close to the radial velocity of the ob- 
ject. Indeed, the observed emission component is relatively nar- 
row and peaks at a velocity very close to that of the SiO maser 
(Fig.EJ. 

This is the first reliable and precise measurement of the radial 
velocity of V838 Mon obtained from optical spectroscopy. The 
result agrees very well with that derived from the SiO maser 
observations. In conclusion, the heliocentric radial velocity of 
V838 Mon is Vh = 71 km s"'. 
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Fig. 5. Sample of the TiO band-heads illustrating the observed 
velocity components marked with dashed vertical lines. The ra- 
dial velocity is in the heliocentric frame. 



3.6. Molecular bands 

The molecular absorption features in the prese nted spectrum ar e 
very similar to those observed in 2005 (.Kaminski et al.ll2009l) . 
We identified bands of TiO, VO, ScO, YO, and AlO. Most of 
the molecular bands are significantly deeper than in standard 
stars of a similar spectral type to V838 Mon. However, in the 
present spectrum, they are mostly ob served to be w eaker than in 
the 2005 one. Similarly as in Kaminski et a n d2009h . when inter- 
preting the observed band shapes, simulations of a homogeneous 
absorbing layer situated above the stellar photosphere were per- 
formed. Usually, a velocity dispersion of 10-20 km s"' has to 
be adopted in the simulations to reproduce the shapes of indi- 
vidual absorption features. Therefore, estimates of the radial ve- 
locity from the bands have a typical uncertainity of ~ 10 km s"' . 
We do not present detailed figures comparing the results of the 
simulations wit h the observed profiles , since they are similar to 
those shown in iKamiriski et all (l2009h . Detailes of the molecu- 
lar bands, for which we were able to derive radial velocities, are 
given in AppendixlAl 



3.6.1. TiO 

As can be seen from Fig. |5] different TiO bands show different 
radial velocities. From the bands with relatively sharp heads we 
were able to identify four velocity components, i.e. -125, -88, 
8, and 88 km s which are marked in Fig. |5] A general trend 
is very similar to that observed in Tvlenda et al. (2009, see their 
Figs. 3 and 4), namely that bands from low excitation levels are 
seen at low (negative) velocities, while the absorptions from ex- 
cited vibrational levels show high (positive) velocities. The only 
important diff'erence is that the highest radial velocity is now ob- 
served to be ~88 km s"' , instead of ~58 km s"' as in 2005. 

The most negative componen t, i.e. -125 km s ', is the same 
as in 2005 (iKamiriski et alJl2009h and is observed only in bands 



Fig. 6. Part of the 2009 spectrum showing details of highly ex- 
cited rotational lines of the TiO y (5,3) and (6,4) bands. The 
V838 Mon spectrum (in black) is overlaid with an M5 III spec- 
trum (in red). The latter spectrum was shifted to a heliocentric 
radial velocity of 88 km s A narrow Ha emission line and a 
broad Ca I absorption profile can be recognized in the spectrum 
of V838Mon. 



absorbing from the ground vibrational state of the TiO system. 
The -88 km s"' component is seen in the bands absorbing from 
the ground vibrational level of TiO, as well as of VO C-X (0,0) 
and ScO A-X (0,0). The component at 8 km s ' is particularly 
well seen in the AlO bands and is discussed below. 

The most redshifted component at 88 km s ' is observed in 
the excited bands, not only in their heads, but also in their re- 
solved rotational lines. An example is shown in Fig. |6] where a 
portion of the V838 Mon sp ectrum is over laid with a standard 
M5 III spectrum (taken from'Bag nulo et aD '2003). An excellent 
fit of the positions of the TiO rotational lines in both spectra was 
obtained when the comparison spectrum was shifted to a helio- 
centric radial velocity of 88 km s ' . This velocity is observed in 
practically all highly-excited (from the higher vibrational levels 
or high-J levels) and weak absorption bands of TiO. Note that 
this is the same velocity as that of the infalling components ob- 
served in numerous atomic lines, as discussed in Sect. 13.31 

Similarly as in Kamiiiski et al ] (l2009h . an estimate of the ex- 
citation temperature of molecular bands can be obtained from 
fitting the results of our simulations to the observed band shapes. 
The shape of the TiO y (1,0) band can be well reproduced by a 
layer of gas excited to ~500 K. To reconstruct the shape of the 
band (2,1) absorbing from the first excited vibrational levels a 
higher value must be assumed, i.e. above 750 K. Thus the exci- 
tatio n conditions of the TiO bands are similar to those in 2005 
(see Kaminski et al. 2009), indicating that the bands showing a 
high blueshift are formed well above the photosphere. A column 
density estimated from the TiO y (1,0) band is now a factor of 2 
lower than in 2005. 
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3.6.2. ScO 

The shape of two components of the ScO A'^II-X'^S^ (0,0) band 
is very similar to that observed in 2005. The simulations show 
that one cannot derive a single value of the excitation temper- 
ature from the observed contour The head formed by lines ab- 
sorbing from highly excited levels, i.e. ^Ric at 6064.31 A, needs 
an excitation temperature above 750 K, and is formed at a ve- 
locity of 66 km s"' (45 km s"' in October 2005). The dominant 
round absorption contours at 6042 and 6082 A, formed by the Q 
and P branches, indicate velocities as low as -10 km s"' and an 
excitation temperature below 350 K. The features of the A - X 
(0,0) band, formed by the '^Ric+'^Qic and ^Qic+^Ric heads, 
have velocities of about -55 km s '. This velocity component 
was also observed in the 2005 spectrum, but now it indicates 
a lower column density. An additional velocity component at - 
30 km s ' can be recognized in the discussed heads in the present 
spectrum. 

The ScO A - X (0,1) band, observed in emission at 6410 and 
6460 A in 2005, is not seen in the present spectrum, neither in 
emission nor in absorption. 



3.6.3. YO 

The YO subband of A^H-X^I.^ (0,0) at 5975 A is well repro- 
duced when assuming an excitation temperature in the range of 
350 - 500 K. The velocity of the absorbing gas is then about 
16 km s '. Similarly, from the fit to the second subband at 
6135 A one can infer an excitation temperature of ~350 K and 
a velocity of 36 km s ' . Both the excitation temperature and the 
colum n density are sim ilar to those obtained from the 2005 spec- 
trum (iKamiriski et alJ l2009i) . However, in 2005, the YO band 
components were observed at a velocity of -8 km s ' . 



3.6.4. VO 

In the observed spectral range only vibrational bands of the VO 
C^^S" - X'^S" electronic system are observed. In two bands, i.e. 
(0,0) at 5735 A and (1,0) at 5468 A, absorbing from the ground 
level, three velocity components can be inferred from their 
shapes. The most blueshifted narrow feature is at -74 km s"'. 
A less clear and dispersed in velocity component is observed in 
a range of -20 to -40 km s"'. Another component dominating 
the redward part of the band has a velocity of 20 km s"' . A less 
intense VO C-X (2,0) head is observed at a velocity of about 
6 km s"', but here also an extension of absorption towards neg- 
ative velocities is evident. 

In 2005, a velocity component of -77 km s"' was observed in 
the (1,0) and (0,0) bands (Kaminski et al. 2009). Within the Hm- 
its of uncertainity we can identify it with the -74 km s"' compo- 
nent in the present spectrum. From our simulations of the (0,0) 
band in 2009, we can infer an excitation temperature of ~300 K, 
i.e. a similar value as that obtained in 2005. However, the column 
density in 2009 decreased by at least a factor of 2 compared to 
the result from the 2005 spectrum. 



3.6.5. AlO 

Contrary to TiO, the observed AlO bands have generally very 
sharp heads if they are absorbing from excited vibrational levels. 
A sample of the AlO band heads is shown in Fig.|7] Our analysis 
of them can be summarized as follows. 




-200 



Velocity (km 



Fig. 7. Detailed view of the AlO B^S - X^E band-heads in the 
velocity scale (heliocentric frame). The characteristic velocity 
of 8 km s ' is marked with the vertical dashed line. The Ha 
emission is also shown at the top of the figure. 



In the spectral range of the present spectrum there is only one 
band of AlO B^S-X^S originating from the ground vibra- 
tional level, (0,0), at 4842 A. The head of this band is blended 
with TiO band rotational lines. Apparently, the band-head 
extends in velocity from the minimum of absorption at 16 
down to at least -100 km s"' . 

All AlO B-X bands originating from excited vibrational lev- 
els have sharp band-heads at ~8 km s with a tendency to 
shift up to ~25 km s ' with the increasing vibrational num- 
ber of the lower level. 

A regular pattern of the rotational transitions formed by 
blends of two main P-branches can be easily recognized in 
the band contour. A good example is the AlO (0,1) band at 
5079 A. The radial velocity obtained from fitting the syn- 
thetized spectrum amounts to 52 km s"'. With the increas- 
ing rotational number J a new component with a velocity of 
~83 km s ' appears, overlapping the blueshifted lines. This 
pattern repeats in the absorptions from the v" =2 level. Above 
v"=2 the rotational structure is consistent with a velocity of 
88kms-'. 

Redward from the band-head, the rotational pattern seems 
to be smoothed out, either by a velocity gradient, or by a 
high density of transitions of the two R-branches, or by both 
mechanisms. 



The rotational contours of the AlO bands are relatively well 
reproduced if the excitation temperature is of ~1500K. A simi- 
lar value was obtained from the AlO bands in the 2005 spectrum 
(Kamiris ki etani2009l) . The most significant difference between 
2005 and 2009 is the position of the heads of the excited vibra- 
tional bands. The present value of 8 km s"' is to be compared 
with -82 km S-' in 2005. 
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3.6.6. Summary 

The main conclusions that we drawn from the above description 
of the molecular bands are: 

- The excitation temperature of the molecular bands of all ob- 
served molecules in 2009 is consistent within the uncertain- 
ties with the values derived from the 2005 spectrum. 

- Differences in the shape of the molecular bands observed in 
the two epochs result mainly from the lower column density 
in 2009, typically by a factor of 2. 

- The kinematic structure as inferred from the heads of the 
molecular bands is generally similar in both epochs. 

- Exceptional are the highly excited lines, which are formed at 
88 km s"' in 2009, compared to 57 km s"' observed in 2005. 

- The change in the velocity is most spectacular in the excited 
AlO bands. In 2009, these heads are at 8kms whereas in 
2005 the head velocities were -89 km s ' . 

4. Conclusions and discussion 

Globally, V838 Mon did not significantly evolve between 2005 
and 2009. It is now slightly more luminous and hotter than it was 
in 2005 (Sect. l3.2T ). It continues to lose matter although at a lower 
rate. This is indicated by the disappearance of clear emission 
components in the atomic lines showing P-Cyg profiles in 2005 
(Sect.[3]l. 

Several atomic absorption lines display a narrow and faint 
emission component. Our analysis shows that this is a real emis- 
sion, most probably a remnant of the P-Cyg emission compo- 
nent observed in 2005, and that its position measures the ra- 
dial velocity of the object (Sect. 13.51 ). The peak of the feature 
is well defined and in all lines, where it is clearly seen, it gives 
Vh = 71kms"' (equivalent to Vlsr - 54kms"'),the sa me value 
as tha t of the SiO masers (Deguchiet al. 2005; Clausse Ti et al.l 
I2007h . A very similar radial velocity was also measured for a 
molecular regio n seen in a near vicini ty of V838 Mon in the CO 
rotational lines (iKaminski et al.ll2bl ih . The region is likely to be 
associated to the dusty matter responsible for the light echo of 
the 2002 eruption. 

An absorption component at radial velocities positive in re- 
spect to the systematic velocity of the object (Sect. 13.31 and 
Fig- [3 is one of the most common features in the atomic lines. 
This component was also observed in 2005 but was less evi- 
dent and seen in fewer lines. A likely reason could have been the 
strong P-Cyg emission component in 2005, which could have 
filled up the redward absorption in certain lines. But it can also 
be argued that the redward absorption component became intrin- 
sically stronger in 2009. This is supported by the appearance of 
a component at 88 km s ' in the highly excited bands of TiO and 
other molecules in 2009 (Sect. 13.611. 

Similarly as in 'Tv lenda et al.l (l2009l) . we interpret the red- 
ward absorption component as evidence of infall motions in the 
atmosphere of V838 Mon. In the present spectrum the maximum 
of absorption is typically observed at Vh ^ 87 km s"', i.e. at an 
infall velocity of ~16 km s The long wavelength wing of the 
feature extends to Vh - 120 km s"', indicating a maximum in- 
fall velocity of ~50 km s"' . This is consistent with the free-fall 
velocity for an 8 Mq star jTvlenda et al.ll2005l) at a photospheric 
radius of ~ 400 Rq (Sect. O, which is -90 km s The in- 
fall motions apparently became more intense in 2009. Note that 
the observed infalls cannot be interpreted as a shrinkage of the 
V838 Mon remnant. With a velocity of 16 km s"' it would result 
in a collapse of the effective radius by a factor of 2 in the course 
of three months. 



The observed kinematics of the AlO bands from the excited 
levels (sharp heads at Vh ^ 8 km s"' - see Sect. 13.6.51 ) is not 
seen in the other molecular bands nor in the atomic line pro- 
files. This suggests that the distribution of aluminium oxide in 
the photosphere and/or the wind is different from the distribu- 
tions of the other molecules and atoms. We calculated distri- 
butions of Al-bearing molecules in a model stellar atmosphere 
taken from the MARCS database (Gustafsson et al. 2008) of an 
effective temperature of 3200 K and log g - 0.0. In the computa- 
tion of the equilibrium composition, equilibrium constants from 
Tsuii ( 1973) were used, except AlO an d AIO2, for which the data 
were taken from a parametrization of ISharp & Huebnerl(ll990l) 
of the JANAF data (Chase 1982). The results reveal that in the 
outer layers of the atmosphere, Al is mainly in the atomic phase. 
AlO is more than two orders of magnitude less abundant than the 
atomic Al, and about two orders of magnitude more abundant 
than AlOH and A1H0 A possible explanation of the observed 
AlO structures is that AlO is additionally produced and exited 
in processes of non-equlibrium chemistry in a postshock region. 
This would explain both the presence of the sharp heads and the 
absorptions from the excited vibrational levels, if the postshock 
temperature and density were sufficiently high. The position of 
Ha emission, as seen in Fig.|7] with its blue edge at ~8 kms"' 
could also have been explained in this hypothesis. 

Absorption features showing the highest outflow velocities, 
i.e. Vh between -130 and -50 km s"', are observed in the molec- 
ular bands arising from the ground vibrational level (Sect. 13.61 
and Table lA.lb and in some atomic lines, mostly in the reso- 
nance ones (Fig. [3J. With the systematic velocity of the object 
(Vh = 71 km s"'), these features imply outflow velocities of 
120-200 km s"'. The structures of these high outflow-velocity 
molecular bands usually give very low excitation temperatures, 
i.e. of 300-500 K (Sect. [fa. 

In 2005, a relatively narrow absorption component (NAC) at 
Vh - -82 km S-' was observed in RbU7800 and MgU4571 
(see Fig. 6 in Kamiriski et al..2009 ). These two lines are, unfortu- 
nately, outside the spectral range of the 2009 spectrum. However, 
a similar component at practically the same velocity is now ob- 
served in Ca I /1 6572 and Li 1/16708 (see panel a in Fig. [3] and 
Fig.U^ 

In lTvlenda et al.l (l2009h . all features in the atomic lines were 
intepreted in terms of a standard spherically symmetric wind 
model. Within this model, the NAC would have originated from 
a local density maximum in the wind, produced by a short-lived 
high mass-loss event. A jet-like structure in the ongoing mass 
outflow was considered as another possible explanation of the 
NAC. A third possibility was that the NAC had nothing to do 
with the current mass outflow from V838 Mon, but instead it was 
formed in a distant matter ejected by V838 Mon during the 2002 
eruption. This last explanation was rejected mainly on an argu- 
ment raised from a discussion of the Til lines. That argument 
was incorrect, however, since the excitation energy of the upper 
level of the T i I line s was mistaken as that of the lower level in 
iTvlenda et al.l(l2009 l). Below we show that the 2002 ejecta is the 
only reasonable explanation not only of the NAC, but probably 
of all the absorption features observed at high outflow velocities. 

A short-lived mass-loss event in the ongoing mass outflow 
can immediately be rejecte d in view of Sect. 13.41 and Fig. |4] In 
the standard wind model of iTvlenda et al.l (120091) . the shell pro- 



^ The difference s with previous computations (e.g. 

iTenenbaum & Ziurvsl l2009h that predicted a much lower concen- 
tration of AlO are caused by the newer equilibrium constants used in 
our modelling. 
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ducing the NAC was at ~2 radii of the central star Therefore, 
the feature should have evolved on a time scale of ~2Q days. 
The features at Vh ~ -50 km s ' in the Ca I and Li I lines were 
very stable over a time span of almost three months. A jet-like 
outflow could in principle produce a stable feature at the pro- 
jected velocity of the jet. However, it is not clear how it could 
produce such a complex and stable structure as that observed 
in the Li I /1 6708 (condensations in the jet would flow out on a 
time scale of days). Besides, as noted above, the molecular bands 
displaying high outflow velocities show very low excitation tem- 
peratures, so they cannot be formed close to the photosphere of 
V838 Mon. 

During the time span between our individual observations in 
2009, matter flowing with a velocity of ~150 km s ' covered a 
distance of -1500 Rq, i.e. ~4 radii of the V838 Mon remnant 
(see Sect. 13.21 ). Therefore, the only reasonable way of explain- 
ing the observed stability of the high outflow velocity features 
in the Ca I /1 6572 and Li I /1 6708, as well as of the narrow ab- 
sorptions at Vh = -92 km s ' in the Bal, Ball, and Til lines 
(see panels b and c in Fig. [3]), is to conclude that the above es- 
timated distance of ~1500 Rq is negligible compared to the dis- 
tance of matter in which the features are formed. On the other 
hand, comparing the line profiles observed in 2009 with those 
observed in 2005, we see that the discussed features evolve on 
a time scale of years. The conclusion is quite clear: the high- 
velocity features are formed in the matter ejected by V838 Mon 
in 2002. During the 2002 eruption the observe d outflow ve- 
locities were between ~ 100 and ~ 600 km s~' dMunari et al.l 
l2002al:ICrause et al.ll200l iKipper et al.ll2004 . The most intense 
mass loss probably took place at the end of the eruption, i.e. 
in March 2002 (Tvlenda 2 005.), with outflow veloci ties ranging 
from -100 to -300 km s"' lCrause et al.ll2003l; iTvle nda 2005|). 
In the 2005 and 2009 spectra, we thus observe relatively slow 
but dense parts of the 2002 ejecta. Faster and presumably less 
dense regions quickly became to o thin to be observable. 

The above hypothesis follows lLvnch et al](l2004r2007h . who 
concluded from their earlier observations that strong molecu- 
lar bands observed in the near-IR were formed in the expand- 
ing 2002 ejecta. It explains several observational facts from our 
observations that are not easily understood otherwise, i.e. abnor- 
maly strong molecular bands in the spectrum of V838 Mon when 
compared to typical giants or supergiants of a similar spectral 
type, very low excitation temperatures of the molecular bands 
formed at high outflow velocities, and decrease in the column 
densities observed in the high velocity molecular bands between 
2005 and 2009 (Sect. l3.6T ). It is now clear also why the maximum 
observed outflow velocity of 215 km s"' was the same in 2005 
and 2009 and removes problems with explaining that high value 
in terms of the winds from cool giants (Tvlenda et al. 2009). 

From the l ow-excitation TiO and VO bands observed in the 
October 2005. lKamiriski et alj (l2009h derived a column density 
of H atoms of ~24 dex cm~^, i.e. ~2 gcm~^. If the matter had 
been lost in March 2002 with a velcity of 150 km s"', it would 
have been at a distance of ~ 2.5 x 10"* Rq in October 2005. 
Assuming spherical symmetry, one obtains a mass of the ex- 
panding shell of ~ 0.04 Mq. This result agrees well with an es- 
timate of the expanding mass made in Lynch et al. ( 2004) and 
is within the limits of the mass los t by V838 Mon in the 2002 
eruption obtained in lTvlendal (l2005l) . 

The disappearance of the B3 V companion from the spec- 
trum of V838 Mon is the most intriguing event in the recent 
years' history of the object. It start ed with an eclipse-like event in 
November-December 2006 (Mu nari et al.l i2007). Then the com- 
panion reappeared in the spectrum for a few months and next 



started a slower but deeper and more complex decline. Note that 
after the 2002 eruption the U brightness of V838 Mon was dom- 
inated by the B3 V companion. Our spectroscopy shows that 
in 2009 the B3 V companion was at least 30 times fainter in 
th e optical than i n 2005 (Sect. lTTT ). Following the interpretation 
of ITvlenda etal] (120091) . we can conclude that the B3 V com- 
panion is now engulfed in a dense dusty cloud of matter lost 
by V838 Mon during the 2002 eruption. The luminosity of the 
companion is thus expected to be now reemitted in the infrared. 
Therefore, one should look for traces of the companion in IR ob- 
servations. This can appear to be a hard task, however, since the 
luminosity of the B3 V companion was only of a few per cent of 
that of the V838 Mon remnant. 

Given that the B3 V companion has already disappeared for a 
few years, one can speculate that in 2007 the companion entered 
the V838 Mon remnant and is now completely engulfed in the 
remnant itself. This idea would imply a second merger phase in 
the history of V838 Mon, in which the B3 V companion would 
finally merge with the core of V838 Mon. This would result in 
a much more energetic and dramatic event than the 2002 erup- 
tion. We do not follow this speculation any further, as it does not 
seem very likely in view of the observational facts. Since its dis- 
covery in October 2002 until the eclipse-like event in 2006, the 
spectrum of the B3 V companion showed no interactions of the 
star with the 2002 ejecta. In particular, the U brightness, domi- 
nated by the B3 V companion in that epoch, remained virtually 
constant. ' This shows that the companion is at a significant dis- 
tance from V838 Mon, so that the matter ejected in 2002 was 
not able to interact with the companion until it reached the star 
at some point in 2006. 

Future observations will distinguish between the two scenar- 
ios. It is very likely that V838 Mon has not yet finished to suprise 
us. Undoubtedly, the object deserves further observational mon- 
itoring. 
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Appendix A: Molecular bands in the spectrum of 
V838 Mon observed in 2009 

A detailed ide ntification of th e molecular features in the spectrum of V838 Mon 
was made in iKamihski et a l. 2009, see their Table 3). The list of the molecu- 
lar bands observed in 2009 remains practically the same as in 2005. Table IaTT] 
presents details of those molecular bands observed in the spectram of V838 Mon 
in 2009 for which we were able to estimate the radial velocity. First column dis- 
plays the observed wavelengths of the bands, while their laboratory wavelengths 
are given in the second column. Identification details of the bands are presented 
in the third column. Heliocentric radial velocities of the bands are given in the 
fourth column. Three different methods were used to derive the velocities, as 
indicated in the last column of the table. The meanings of the symbols are as 
follows: M5 in - comparison of the V838 Mon spectrum to the UVES spectram 
of an M5 III star HDl 18767 i Bagnulo et al. 2003); LAB - direct comparison of 
the band head position with the laboratory measurements (see Kaminski et aP 
[2009, for la boratory data sources); SYNTH - spectral synthesis of the bands as 
described in lKamihski et aU j2009l) . 
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Table A.l. Heliocentric radial velocities of the molecular bands in the 2009 spectrum of V838 Mon 



'Jobs 


'llab 


Identification 


Velocity 


Method 


(A) 


(A) 




(kms-') 




4788-4804 




TiO a (2,0) high J tail 


+88±5 


M5 III 


4805.2 


4804.333 


TiO a (3,1) R2 


+56±25 


LAB 


4831-4841 




TiO a (3,1) high J tail 


+88±5 


M5 111 


4842.47 


4842.27 


AlO B^S-X^S (0,0) R2 


+12±10 


LAB 


4889.26 


4889.0 


AlO B^S-X^S (2,2) R2 


+16±10 


LAB 


4953.0±0.2 


4954.592 


TiO a (1,0) 


-90±18 


LAB 


4993-4998 




TiO (1,0) high J tail 


+88±10 


M5 III 


5042-5078 




TiO a (3,2) 


-i-88±10 


M5ni 


5079.5 


5079.36 


AlO B^X-X^E (0,1) 


-i-8±5 


LAB 


5087-5096 




AlO B^S-X^S (0,1) P-branch 


-i-52±5 


SYNTH 


5102.3 


5102.13 


AlO B^Z-X^S (1,2) R-head 


-i-10±5 


LAB 


5113-5118 




AlO B^S-X^S (1,2) P-branch 


-i-52±5 


SYNTH 


5123.6 


5123.33 


AlO B^S-X^S (2,3) R-head 


-i-15±6 


LAB 


5143.3 


5142.89 


AlO B^S-X^S (3,4) R-head 


-i-23±10 


LAB 


5130-5137 




AlO B^S-X^S (3,4) high-J 


-i-88±10 


SYNTH 


5164.4 


5166.450 


TiO a (0,0) R2 


-119±10 


LAB 


5197-5203 




TiO a (0,0) high-J 


-i-88±10 


LAB 


5213-5225 




TiO a (0,0) high-J 


-i-88±10 


LAB 


5226.3 


5228.2 


VO C-X (2,0) R-head 


-109±10 


LAB 


5228.0 


5228.2 


VO C-X (2,0) R-head 


+11±10 


LAB 




5240.463 


TiOf'A-a' (0,0) 


+88±10 


M5 III 




5239.418 


TiO a (2,2) R2 


+88±10 


M5 111 




5275.8 


VO C-X (3,1) 


+88±10 


M5 III 




5307.116 


TiO a (3,3) 


-l-88±10 


M5 III 


5377.21 


5376.813 


AlO B^Z-X^S (2,4) 


-i-22±10 


LAB 


5467.7 


5469.3 


VO C-X (1,0) 


-88±10 


LAB 


5468.3 


5469.3 


VO C-X (1,0) 


-52±10 


LAB 


5488-5498 




TiO a (0,1) high-J 


+88±10 


M5 III 




5496.423 


TiO a (1,2) R2 


+88±10 


M5 III 


5512-5525 




TiO (1,2) high-J 


+88±10 


M5 III 






TiO fi (0,0) 


+88±10 


M5 III 






TiO /3 (1,1) 


+88±10 


M5 III 






TiO ji (2,2) 


+88±10 


M5 III 


5650-5660 




TiO 7 (3.1) high-J tail 


+88±10 


M5 III 




75712.6 


TiO Y (4,2) 


+88±10 


M5 111 




75717.0 


TiO y' (4,2) 


+88±10 


M5 III 


5735.22 


5736.703 


VO C-X (0,0) 


-78±10 


LAB 


5735.98 


5736.703 


VO C-X (0,0) 


-38±10 


LAB 




5758.741 


TiO a (0,2) R2 


+88±10 


M5 111 


5845.73 


5847.593 


TiOr' (1,0) Ri 


-97±10 


LAB 


5950 




TiO 7' (2,1) R3 


+88±10 


M5 III 


5999.64 


6001.151 


TiOy (3,0) R2 


-75±10 


LAB 




6064.31 


ScO (0,0) Ri 


-i-88±10 


M5III 


6077.65 


6079.30 


ScO (0,0) Qi 


-81±20 


LAB 


6147.15 


6148.68 


TiO Y (0,0) ■^■R2i 


-75±10 


LAB 


6269.3 


6268.86 


TiO Y (1,1) 


+18±10 


M5 III 


6292.2 


6294.80 


TiO y (2,0) R3 


-123±20" 


LAB 


6294.2 


6294.80 


TiO y (2,0) R3 


-29±10'' 


LAB 


6296.0 


6294.80 


TiO y (2,0) R3 


+57±5'' 


LAB 


6318.9 


6321.21 


TiO y (2,0) R2 


-110±10'' 


LAB 


6320.7 


6321.21 


TiO y (2,0) R2 


-24±10° 


LAB 


6349.7 


6351.29 


TiO y (2,0) Ri 


-75±20" 


LAB 


6398-6412 




TiO y (3,1) high-J tail 


+88±10 


M5 III 




6447.90 


TiO y (4,2) 


+88±10 


M5 III 




6478.55 


TiO y (4,2) Ri 


+88±10 


M5 III 




6512.409 


TiO y (5,3) R2 


-i-88±10 


M5in 




6543.950 


TiO y (5,3) Ri 


-i-88±10 


M5in 




6549.600 


TiO y (6,4) R3 


-i-88±10 


M5in 




6577.716 


TiO y (6,4) R2 


-i-88±10 


M5in 


6649.34 


6651.252 


TiOy (1,0)R3 


-86±20 


LAB 


6678.81 


6680.796 


TiOy (1,0)R2 


-89±20 


LAB 


6712.59 


6714.477 


TiOy (1,0) Ri 


-84±20 


LAB 




6717.578 


TiOy (2,1) R3 


-h68±5 


M5m 


6733-6740 




TiOy (1,0) high-J tail 


-h68±5 


M5m 




6747.589 


TiOy(2,l)R2 


-h68±5 


M5in 


6765-6780 




TiOy (2,1) high-J tail 


-h88±5 


M5in 


6808-6815 




TiOy (2,1) high-J tail 


-h88±5 


M5in 



"No clear bandhead, a velocity gradient in the profile is observed. 



